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In J. Math. Phys. 13, 1608-1621 (1972), Erdahl considered the convex structure 
of the set of iV-representable 2-body reduced density matrices in the case of fermions. 
Some of these results have a straightforward extension to the m-body setting and 
to the more general quantum marginal problem. We describe these extensions, but 
can not resolve a problem in the proof of Erdahl's claim that every extreme point is 
exposed in finite dimensions. Nevertheless, we can show that when 2m > N every 
extreme point of the set of iV-representable m-body reduced density matrices has a 
unique pre-image in both the symmetric and anti-symmetric setting. Moreover, this 
extends to the quantum marginal setting for a pair of complementary m-body and 
(N — m)-body reduced density matrices. 
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I. INTRODUCTION 

The development of quantum information theory has generated interest in what is known 
as the quantum marginal problem, i.e., the question of when a given set of reduced density 
matrices (RDMs) can be obtained from an iV-body state. A special case, known as the N- 
representability problem asks when an m-body RDM for a fermionic system can be obtained 



2 



from an iV-body fermion state. 

This question was extensively studied in the 1960's and 1970's in the hope of finding a way 
to use the 2-body RDMs to simplify complex iV-body computations. The iV-representability 
of the 1-body RDM has a simple solution found independently by several groups, including 
[4, 7, 11, 21]. A 1-body RDM is iV-representable if and only if its eigenvalues satisfy the 
condition < < which expresses the Pauli exclusion principle. 

However, iV-represent ability for the 2-body RDM is a far more challenging problem 
on which little progress was made for over 30 years. In 2007, the hardness of this prob- 
lem was formally recognized by proving that it belongs to the complexity class known as 
QMA-complete [15], which is the quantum analogue of NP-complete, i.e. testing for N- 
representability would require exponential time even on a quantum computer in the worst 
cases. 

Slightly earlier, Klyachko [1, 10] found a complete solution of the pure state N- 
representability problem for the 1-body RDM. The earlier solution of the mixed state prob- 
lem can be restated as the fact that the extreme points of the convex set of 1-body RDMs 
arise from iV-body Slater determinants. However, little was known about the images of more 
complex states beyond an abstract induction result of Coleman [4] which can be viewed as 
a constrained version of Weyl's problem [19]. 

These advances came as a result of recent work in quantum information theory in which a 
RDM is often called a "quantum marginal" . There are several variants [3, 9] of the question of 
whether of a given set of quantum marginals is consistent; the iV-representability problem 
can be regarded as a special case of one of these. Although consideration of symmetry 
constraints is at the heart of iV-representability, some fundamental aspects of the convex 
structure carry over from the anti-symmetric iV-representability situation to the more general 
quantum marginal setting for m-body RDMs. 

In [6, Section 3], Erdahl claims that in finite dimensions, every extreme point of the convex 
set of iV-representable reduced density matrices (RDMs) is exposed. In trying to extend this 
argument to quantum marginals, we found a subtle problem in the proof, which we describe 
in Section III. Because some of the results in [6] easily extend to quantum marginals, we 
introduce notation in Section II which makes such generalizations transparent. 

However, Erdahl's primary application of this result, i.e., his proof [6, Section 6] that 
when 2m > N, the pre-image of an exposed point of the set of iV-representable fermionic 
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m-body RDMs is unique holds for arbitrary extreme points. One needs only the observation 
that the pre-image of a face is a face. 

In addition we use one of his key observations [6, lemma 6.1] to show that a pair of 
complementary quantum marginals has a unique pre-image when both are extreme points 
of their respective sets of RDMs. This is a rather remarkable property of extreme points 
because, in general, a pair of density matrices whose Hilbert spaces have a non-trivial overlap 
seems to be essential to the unique determination of a pure state pre-image. Moreover, this 
contrasts sharply with the situation when 2m < N for which a single extreme m-body RDM 
can have multiple pre- images as shown recently in [16]. 

The rest of this paper is organized as follows. Section II introduces necessary concepts 
and notation for both anti-symmetric case and quantum marginals. We then describe the 
problem in Erdahl's proof of the claim that every extreme point is exposed in finite dimen- 
sions in Section III. Inspired by Erdahl's primary application in [6, Section 6] that when 
2m > N every exposed point of the set of iV-representable RDM has a unique pre-image, 
we study the pre-images of extreme points and interior points in Section IV. Specifically, 
before describing our main results of Section IV, we first give some intuition and background 
in Sub-Section IV A. Then, by applying a technical lemma provided in Sub-Section IV B, 
we show that Erdahl's result holds for extreme points even if they are not exposed in Sub- 
Section IV C. In Sub-Section IV D, we show that if the contraction map has a non-trivial 
kernel, then every interior point of the set of iV-representable RDM has non-unique pre- 
images. Conclusion and future work can be found in Section V. 



II. NOTATION 



A. The antisymmetric case 

Let % be a <i-dimensional Hilbert space and let = %® N = T-L AN be the anti-symmetric 
subspace of its iV-fold tensor product. Let T>(X) denote the set of density matrices for an 
antisymmetric subspace X C %® N , i.e., 

V(X) = {p E B(X) : p > 0, Trp=l} (1) 

where B{X) is the set of all bounded linear operators on X. 



i 



Let T> m (X) denote the corresponding set of reduced density matrix (RDM), i.e., 

V m (X) = {p m G V(U® m ) : 3 p G V(X) with Tr m+1 _ N p = p m } 

= convex hull {Tr m+1 ,... )JV |V>)(V>| : \if>) E X}. (2) 

Then P m ("H AAr ) is the set of iV-representable RDM. The polar cone is 

V[V m (X)] = {V e B(H® m ) : Tr (Vp m ) > V p m G £> m (X)}. (3) 

Every V G V[V m (X)] can be associated with a positive semi-definite Hamiltonian in B(X) 
given by 

*W0= £ V{jlj2...,jm) (4) 

Jl.J2,-jm 

where V(ji,j 2 ■ ■ ■ ,j m ) acts on ® Hj 2 <g) . . . ® "H^. An extreme point is exposed if there 
is a 1/ G P[X? m (^ AAr )] such that Tr{Vp m ) = and Tr (V-? m ) > for all 7 m ^ p m . It 
then follows immediately from the variational principle that the pre-image of p m is a density 
matrix for the ground state eigenspace of Hn(V). (Although the generic situation is a one- 
dimensional eigenspace, it was shown in [16] that the pre-image of an exposed point can 
be the set of density matrices associated with a degenerate ground state eigenspace of an 
m- local Hamiltonian.) 

Let {p\ = Xp + (1 — X)p' : A G [0, 1]} C C denote a line segment of a convex set C. A face 
J 7 of C is a subset with the following property. Whenever some element p\ of a line segment 
is in T for some A G (0, 1), then the entire line segment is in J 7 . An extreme point of C is 
a face of dimension 0. A face J 7 of T> m (X) is exposed if there is an element V of the polar 
cone such that Tr (Vp m ) = Vp m G T and Tr (V^ m ) > for all 7 m ^ J 7 . 

Note that if a density matrix p G £*(A) does not have full rank, then it is in the face 
(ker / o)" L nP(X), i.e., the set of all density matrices 7 for which range 7 C range p. Therefore, 
any state in the interior of T>(X) must have full rank. 

B. Quantum marginals 

Some results about the convex structure of the set of ^-represent able RDMs have 
straightforward extensions to convex sets associated with quantum marginal problem. In 



5 



this situation we want to know if there is an iV-body state consistent with a set of m-body 
RDMs. 

In the general case, we can replace T-L AN by an arbitrary subspace Hn Q %® n '. For a 
given set of indices J = {ji, J2, ■ ■ ■ ,jm}, let J c denote the complement in {1,2,..., N} and 
for any iV-body state pi,2,...,N G T^{Hn) define 

Phh-im = Pj = Ttjcpi,2,...,N- (5) 

What we called p m in the previous section is more properly called p\ > 2,..., m - When 7-L N is the 
symmetric or anti-symmetric subspace of 7-l® N this slight abuse of notation is justified by 
the fact that p m determines Pj 1 j 2 ...j m f° r anv set of indices. If we let Hj C denote the 
subspace induced by if at, then, as above for any subspace X C H N we define 

Vj{X) = {pj G V{Uj) : 3 p 1)2) ... !jv G V{X) such that pj = TWi,2,..,iv}- (6) 

We can then replace p m by a vector, e.g., 

R2 = {Pl2, Pl3, ■ ■ ■ 1 Pin, P23, ■ ■ ■ 1 PN-1,n) 

or, more generally 

R-m — (Pl,2,...,mj Pl,2,...,m-l,m+l> • • • 5 PN-m+l,...,N-l,N) 

consisting of all possible m-body RDMs in some prescribed order. Then we define 

V m (X) = {R m G V(X) : 3 p G V(X) with Pj G 2?j(X) VJ with |J| = m}. (7) 

Given a vector i? m with elements in T>^H® m ) the consistency problem for quantum marginals 
asks if there is a state pi,2,...jv £ 2^(%v) whose RDMs are given by the elements of R m . 
When "Hjy = "H A7V this is the iV-representability problem. Both T>j{X) and P m (X) are 
closed convex sets. 

The polar cone 7 3 [X ,m (?^jv)] consists of vectors 

Kn — (^1,2, ...,mi ^1,2,... (8) 

with elements \ }.;.,...,-,. G B(H h ® U h ® . . . <g> H im ) and Tr (V m • i? m ) = Y.\j\= m (J T v JPj)- 
Elements of the polar cone of T>{l-i® m ) are associated with m-local iV-body Hamiltonians 
Hn(V) = ^2\j\ =m Vj which are positive semi-definite; the extreme rays have a ground state 
with eigenvalue zero. 

We may and do carry over the descriptions of face and exposed face using the concept of 
polar cone to the quantum marginal case. 
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III. A SUBTLE PROBLEM IN ERDAHL'S LEMMA 



We now describe the problem in Erdahl's proof of the claim that every extreme point 
is exposed in finite dimensions. For simplicity, we consider the original anti-symmetric 
situation with Hn = "H A7V and use the notation introduced in Section II A. However, it 
should be clear that everything goes through in the more general quantum marginal setting 
with p m replaced by R m . 

Let A,B be subspaces with "H A7V DiDB, and suppose that T = V m (A) is an exposed 
face of T> m {T-L hN ) the set of iV-representable RDM and Q = T> m (B) is another face which 
lies in J 7 . In [6, section 3] 1 Erdahl states the following as Lemma 3.1. 

Statement 1 If Q = T> m (B) is an exposed face of the convex set T> m (A), then it is also an 
exposed face of D m ('H AAr ) . 

By assumption, 

a) there is & V £ V[D m (H AN )] which exposes J 7 , i.e., Tr(Vp m ) = V p m E J= and 
Tr(Vp m )>0 V Pm eV m (A)\T. 

b) there is a If £ V[V m {A)] such that Tr (W p m ) = V p m e Q and Tr (Wp m ) > 
W Pm eV m (A)\g. 

Although Hn(V) is positive semi-definite on H AN , the operator Hn(W) need not be 
positive semi-definite. However, Erdahl claims that one can find a t > such that tH^{V) + 
Hn(W) is positive semi-definite. But we claim that this is true if and only if B C ker H^{W). 
(Erdahl states that "clearly" B C ker Hn(W), but gives no proof.) 

We can decompose 

u = Ui ® n 2 © n 3 = b © (a n b l ) © a ± , 

and write Hn(V) and Hn(W) as block matrices accordingly 

/n n n \ 



H N {V) 






A 33 



/ n n v \ 



H N (W) 



Y 13 
Y22 Y23 
Y31 Y32 V33 



(9) 



1 It is ironic that there is a typesetting error in the title of this section so that what was intended as "In Finite 
Dimensions ..." appears as "INFINITE DIMENSIONS ..." rather than "IN FINITE DIMENSIONS" . 
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The matrices X 33 and Y22 are strictly positive definite by construction. However, the 2x2 
block submatrix of tH^iV) + H^iW) obtained by omitting the second row and column is 



Y 13 

Y3I tXsz + I33 

which can never be positive semi- definite for any t > unless Yi 3 
to B C ker H N (W). Although 



(10) 



0, which is equivalent 



EB Tn,... )m ^(Tr 



■m+l,...,7V 



this does not imply that \ip) G ker.ffjv(W) because H N (W) need not be positive semi- 
definite 2 . We can only conclude that G ker PaHn{W)Pa where Pa is the projection onto 
A = ker Hn{V). In general, we do not expect that PaHn(W)Pa is an m-body Hamiltonian. 

If it were true that Yi 3 = 0, we could simplify Erdahl's argument slightly because it 
suffices to consider the non-zero 2x2 submatrix of 



tH N (V) + H N (W) 





Y22 Y23 

Y32 tX 33 + I33 



\ 



(11) 



As noted above, the assumptions on the regions for which Tr (Vp m ) > and Tr (Wp m ) > 
imply that both X33 and Y22 are positive definite. Therefore, using a standard result for 
2x2 block matrices, (11) is positive semi-definite if and only if 



^23^22 X ^23 < tX 33 + Y33 



(12) 



which always holds for t sufficiently large. To be precise, we could find p, A > such that 
Y22 > pi and X 33 > XI so that it suffices to choose 



/^33||). 



(13) 



2 Indeed, if Hn{W) is not positive semi-definite, then there will always be vectors in T-l AN such that 
{tjj, Hn(W)iP) = but ip £ ker Hn(W). This is because on Hjy(W) will have both positive and negative 
eigenvalues on [ker Hn(W)] ± which implies that in the numerical range of Hn(W)\„ Hn ^ w ^ ± . 
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IV. PURE STATE P RE-IMAGES AND THE UNIQUENESS QUESTION 

A. Introduction and Intuition 

In [6, Section 6] Erdahl proves that when 2m > N every exposed point of the set of N- 
representable RDM has a unique pre-image. We show that this holds for extreme points even 
if they are not exposed by showing that [6, Lemma 6.1] still holds in the form of Theorem 6 
below. Remarkably, we use this result to extend Erdahl's result to complementary pairs of 
quantum marginals, even without permutational symmetry or overlap. 

Theorem 2 If 2m > N any extreme point of the set of N-representable RDM has a unique 
pre-image. Moreover, for fixed J whenever both pj and pjc are extreme points of T>j(T-Ln) 
and VjcIHn) respectively, then the pre-image of the pair (pj,pjc) is unique. 

The proof of Theorem 2 is an immediate consequence of Theorem 6 which is a straight- 
forward generalization of Lemma 6.1 of Erdahl [(>]. Before describing these results, we give 
some intuition and background. 

The intuition underlying this result comes from the so-called Schmidt decomposition in 
the special case that p m has non-degenerate eigenvalues with eigenvectors {|Xfe)}- When the 
pre-image is anti-symmetric and m > N — m, we can find pN-m and its eigenvectors {|0fc)}. 
Then any pure state pre-image must have the form 

k 

If is not unique (up to an overall phase factor), then, after suitable absorption of some 
phase factors into the \Xk), one can find a P arr of pre-images which can be written as non- 
trivial superpositions of the form 

|^-> = a#i) + e^j/l^j) = \Xk ® <f>k) + ^ \i k \\k ® 4>k) 3 = 1,2 (15) 

with K\,Ki non-empty and disjoint, |x| 2 + |y| 2 = 1 and < u\ ^ u 2 < 2ir. Since Hn is a 
vector space, both — \ty 2 ) and e~ luJ1 \^i) — e~ luJ2 are in Hn which implies that both 
and ^2) are in %n- Then using the orthogonality of the we find 

p m = |a;| 2 Tr m+ i r .. ;JV (i>i\ + \y\ 2 Tr m+lj ^ N \if) 2 )(i/}2\- (16) 
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This is a mixture and hence not extreme. This argument is essentially due to D. Smith [20] 
who considered the case m = N — 1 and for m = 2, N = 3 even analyzed the general case 
with degenerate eigenvalues. This argument is also the basis for a result of Diosi [5] which 
we state next. 

Theorem 3 (Diosi-Smith) Almost every pure state G 1-1® N is uniquely determined by a 
pair of RDM (pj, pj,) with J n J' ^ and J U J' = {1, 2, ... N}. 

Sketch of proof : For N > 3, the set of \ip) for which pj has degenerate eigenvalues has 
measure zero. Moreover, the hypotheses of the theorem imply that J c C J' so that we can 
determined pjc by taking a suitable partial trace of pji. Therefore, we can assume that 
we have an expansion as in (14). Since we also have {J') c C J the eigenvectors of pj> and 
p( ji)c give a second expansion. Equating these two expansions pointwise gives a set of linear 
equations 3 for Xj = e i6,J and x'j = e te j . We have at most dim % j + dim Hjc unknowns and 
dim "H^r equations. In typical situations dimT-L N « dim %j dim %jc . When this is not true, 
as in the antisymmetric case, the symmetry restrictions give additional equations. Thus, in 
general, we expect a unique solution for Xj and x' k . With the additional requirement that 
acceptable solutions must satisfy \xj\ = \x' k \ = 1, the situations in which multiple solutions 
exist for the phases Uj will be expected to be very rare. QED 

Theorem 2 says that when both pj and pjc are extreme points of their respective set of 
RDMs, then the pair (pj, p ji) determines a unique pure pre- image \ip) although for J 1 = J c , 
J n J' = precluding a second expansion. In view of the argument above, the claimed 
result might seem too good to be true. However, when Hn = H® N with no constraints, 
Vj(U m ) = V(H n <g> H h <g> . . . <g) %J which is simply the convex hull of projections onto 
pure states |%) G l-i® m . Therefore, in the complete absence of constraints, Theorem 2 
simply reflects the fact that the extreme pairs have the form |x){xl; 10) (01 an d the pre- 
image = |x® <fi) is a pure product state. On the other hand, when Hn is the symmetric 
or anti-symmetric subspace of H® N and 2m = N, we also have J = {1, 2 . . . m} and J' = 
J c = {m + 1, m + 2, . . . ,N}. However, permutations of the single expansion (14) which 
exchange j <->• k with j < m and k > m give linear equations which determine the phases, 
as above. 

3 In general, an arbitrary pair of RDM p,j and pjc , will not even have a consistent set of linear equations 
with a solution for Xj and x'p much less a set with unit norm. 
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More generally, one might have other symmetry constraints on Hn, e.g., rotational sym- 
metry or translational symmetry for spin lattices. In such situations, the equations for 
determining the phases might not be as transparent as for permutational symmetry. The- 
orem 2 says that, nevertheless, the condition of being an extreme point of T>j jc(Hn) is 
sufficiently strong to uniquely determine the pre-image. To prove this, we need some addi- 
tional lemmas. 



B. Key Lemmas 

Lemma 4 Let Q be the span of the extreme points of a face ofDCHx) for which every state 
maps to a unique m-body RDM, i.e. 

pj = Tijc\^)^\ V |V) e g. (17) 

Then for any pair of orthogonal vectors \ipk) Q , 

(i) The (N — m)-th order transition density matrix Trjcl^)^^! = 0; 
(in particular, for symmetric or anti- symmetric Tr m+ i ) ...jv|V'j)(^ ; A;| = ®)- 

(ii) For any m-body operator 5 Jlv .. Jm = Bj acting on 'Hj 1 ®T-Lj 2 £§> . . . ®"Hj m , (ipj, Bjipk) — 

0. 

Proof : When dim(? = 1, the result holds trivially, since there are no orthogonal pairs of 
vectors. Therefore, we assume that dim(? > 2. Let {pj} be the unique set of m-body RDM 
onto which all normalized vectors in Q are mapped. Then for any a, b > with a 2 = b 2 = 1 
and any 9, consider the RDM pj of a\ip\) + e b\ip2)' 

pj = Tijc [(a|^ 1 ) + e^|V 2 ))(a(V 1 |+e-^ 2 |) 

= a 2 pj + e ie abTr jC \if)i)(if) 2 \ + e^abTijc |^ 2 )(Vi| + b 2 Pj 

= pj + ab^Tijcl^iH+e^TTjcl^dj^ (18) 



which implies 



= e ie Tr jC I^X^I +e- i6 Tr: J c \^){^\ (19) 



which is equivalent to e l2d A = A< when A = Tr jc \ipi){ip2\- Since this holds for 6 arbitrary, 
A = Tt jC 1^)^1 = 0. 
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To prove (ii) it suffices to observe that (i) implies 

tyuBjfa) = TrjBjTrja^)^ = TrjBjO = 0. (20) 

Erdahl applied this Lemma to exposed points with Q the ground state eigenspace of an 
m-local Hamiltonian. For arbitrary extreme points the existence of Q follows from the fact 
that the pre-image of a face is always a face. For completeness, we sketch an elementary 
proof for extreme points. 

Lemma 5 Let |^>i), ^2) be a pair of orthogonal vectors in H such that Trjc|^i)(^i| = 
Tr jc I ip2 ) ("02 1 = Pj and assume that pj is an extreme point ofDjil-L^). Then TTjc\ip)(ip\ = 
Pj for any unit vector G span{\ipi), \ip 2 ) }■ 

Proof: Let a, b G C, satisfy \a\ 2 + \b\ 2 = 1 and define \if}±) = a\ipi) ± b\ip2). Then it is easy 
to verify that 

|Tr ja |^+) 1 + ±Tr jC |^_) | = pj (21) 

contradicting the assumption that pj is extreme unless Tr jc\ip + )(i/) + \ = Tr jc |^_)(^_| = pj 
which implies 

abTTjc\ip + )(ip^ \ + abTYjc\ipJ)(ip + \ = 

which implies that \i/)±) G span{|^i), \ip2)}- QED 

Note that since (a, b) is an arbitrary pair of complex numbers with \a\ 2 + \b\ 2 = 1, we can 
further conclude that (as above) Trjc\ip + )(tp_\ = 0. 

Theorem 6 (Erdahl) Let {\ipk)} be an orthonormal basis for Q with common m-body RDM 
pj as in Lemma 4- Then for any pair with j 7^ k, range(Trj|-0j)(^ :; |) C ker(Trj| , 0/ c )('0fc|) 
or, equivalently in the case J = {1,2,..., m}, 

\1>i) = X>IXi®«?) ( 22 ) 
t 

with p t and {\xt)}the eigenvalues and eigenvectors of p m and (6^,6^) = 5jk for all s,t. 

Proof: The expansion (22) is the standard Schmidt decomposition with {]#;?)} the eigenvec- 
tors of Trjc\ijjj){ipj\. There is no loss of generality in assuming that p t > 0. By Lemma 4(i) 

= Tijc\i)j){il) k \ 

= J2J2vs»tTr J c\xs®0i){ X t®9't\ 

s t 

= £Z>.M&#>lx.><x*l- (23) 

s t 
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Since the set {|x s )(Xt|} is an orthonormal basis for BCH® m ) with respect to the Hilbert- 
Schmidt inner product, the coefficients in (23) above must be zero. By assumption, the 
singular values p s p t > are all non-zero. Therefore 

(9i,9t) = V s,t when j^k. (24) 

This gives (22) and implies range (Trj|^-) (ipj\) C (ra.iage(Trj\ipk) (ipk\)) ± ■ The result then 
follows from ker A = (rangeA) -1 for any self-adjoint A. QED 

To see how remarkable Theorem 6 is, consider the case N = 5, m = 2 and suppose that 
\ip2) have the same 2-body RDM. Then (22) becomes 

hM = ^2fJ.t\Xt(l,2)® ^(3,4, 5)) 

|^ 2 ) = ^/it|x i (l,2)®0 t (3,4,5)) 
t 

with {\0 t )} and spanning orthogonal subspaces of "H® 3 . However, 

\th\*T± 3 \e t )(e t \ = P 4 5 = Y1 H 2 Tr 3 |0*>(0 t | (25) 
* t 

so that the convex hulls of the 2-body RDM associated with these orthogonal subspaces are 
identical. 

C. Pre-images of extreme points 

When an extreme point is exposed, it follows that the pre-image is the convex hull 
of unit ball of the ground state eigenspace so that the hypothesis of Theorem 6 holds. 
Lemma 5 implies that Theorem 6 also holds for arbitrary extreme points, enabling us to 
prove Theorem 2. 

Proof of Theorem 2: This is an easy corollary of Theorem 6. Let Qj, Qjc be the subspaces 
of "Htv with common RDM pj and pjc respectively. Q = QjV\ Qjc By assumption Q is not 
empty. If dim(? > 2 then it contains a pair of vectors satisfying Theorem 6. But then 
Pjc = Tr = TtjI^jX^I which implies that 6 J S = 9^ contradicting the strong 
orthogonality condition (24). When Hn is the symmetric or anti-symmetric subspace of 
'H <g ' N , and 2m > N, then the same argument can be applied with J = {1,2,..., m} and 
J c = {m + l,...N}. QED 

In the anti-symmetric case, particle-hole duality gives the following 
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Corollary 7 For any N-fermion system, any extreme point ofV m ('H AN ) has a unique pre- 
image if 2m > d — N . 

In view of Theorem 2, and experience with situations encountered in atomic and molecular 
problems, Erdahl conjectured that the pre-image of every extreme point of T> m {T-L hN ) is 
unique. 

In [16], it was shown that this conjecture is false, even in the case of exposed points, 
for both the original fermionic situation and the quantum marginal problem. As observed 
in [16], Lemma 4 implies that a subspace of that form is a quantum error correcting code 
(QECC) which can correct [^^J errors or detect m — 1 errors. Although most QECC are 
associated with interior rather than extreme points of the m-body RDM, examples of QECC 
which arise as the ground state eigenspace of an m-local Hamiltonian are now known. 

These codes are typically associated with spin systems. A standard second quantization 
argument maps the spin eigenstates to fermionic states with only half-filled shells. To ensure 
that these are the ground states, the Hamiltonian includes a "penalty term" of the form used 
in the Hubbard model. The classic examples of extreme points of 2-body RDM exhibit strong 
pairing, as the case of "tight binding" and unique pre-images. 

There are presumably large families of extreme points intermediate between this cases, 
which are yet to be found. 

D. Pre-images of non-extreme points 

Since dim'H®^ is much larger than dim"H® m , the partial trace operation will, in gen- 
eral, take many distinct states to the same output state. Intuitively, this suggests that 
non-extreme points of T> m {l-iN) to have multiple pre-images. Although this is the generic 
situation, there are some significant exceptions which we discuss before stating a result about 
non-uniqueness. Building on notation introduced previously, (e.g., the pair (pj, pj/ ) in The- 
orem 3) we let Q = (p j i; pj 2 , . . . , pj v ) denote a vector of RDM and T>j lt j 2 ^j v (X) the convex 
set of such vectors with pre-images in X. We are interested in the conditions under which Q 
does or does not have a unique pre-image in T>(T-Ln), and let T : V^H^) >->■ 2^j 1 ,j 2i ...,j, / (?^jv) 
for which T(p) — Q (where we have suppressed the dependence on Jj for simplicity). 

If a RDM vector Q does not lie on a face of I , j 1 .j 2 ....,j^('Hv) which is also simplex, it 
will have multiple decompositions into extreme points which we can write Q = J2k a kQk = 
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Ylik a 'kQ'k an d assume that each Q k or Q' k has a unique pre-image \4>k){ipk\ or \ip' k ){ip k \, 
respectively. Then both p = J^fe °fclV , *;)('0fc| an d p' = Ylik a k\ 1 Pk)( 1 Pk\ are pre-images of Q 
which suggests non-uniqueness. However, it is possible that the two states p = p' despite 
the formally distinct decompositions. An example of such a situation occurs when an m- 
local Hamiltonian has a doubly degenerate ground state eigenspace; then the corresponding 
set of iV-body RDM for the ground state is isomorphic to the Bloch ball. The image 
of the Bloch ball under partial trace is also an ellipsoid which corresponds to a face of 
T> m '(Hn) ■ Each element of this face, even those which are not extreme, has a unique (possibly 
mixed) state pre-image. This example also describes a situation in which the map from 
T : V n (Hn) i—t- T> m {Kj^) is one-to-one. One simply chooses the subspace "Hat to be the 
ground state eigenspace of such an m-local Hamiltonian with degeneracy. 

A less artificial example of a one-to-one mapping occurs when Hn = H AN is anti- 
symmetric and dim'H = d = N + m so that dim"H AAr = (^) = f ) = dimH Am . Then 
r : p N h-> p m is given by particle-hole duality as described in [6, Section 4] and [18]. More- 
over, the maps from p N y p m must also be one-to-one for m > d — N . 

In view of these examples, the following theorem seems to be best possible. For the 
anti-symmetric case, a stronger result was claimed by Rosina [17], who apparently did not 
realize that the kernel of V could be trivial. His argument was similar to that given below. 
Since 

ker(r) = {T G B(U N ) : TrjcT = 0, £ = 1, 2 . . . i/}, (26) 

it follows that TrT = 0. Moreover, a non-trivial kernel contains at least one self-adjoint 
matrix which we can write as T = T + — T_ with both T + and T_ positive semi-definite. For 
any p of full rank and sufficiently small p, the matrix p+pT £ T>{T-Ln) and T(p+pT) = T(p). 
Thus, T has a non-trivial kernel if and only if 2?j 1) j a ,... ) j„('Hjv) contains two distinct density 
matrices with the same pre-image. 

Theorem 8 If the map T : p i— » Q has a non-trivial kernel, then all interior points of 
/ Dj 1} j 2 ,... i j 1/ (Hn) have non-unique pre-images in T>(Hn). Moreover, in this situation, every 
interior point has at least two pre-images on the boundary of 12(1-1 n)- 

Proof: Let Q be an interior point of I'j 1 ,j 2 ,...,j 1 ,(^jv), Since T is linear 4 the inverse image of 



4 Here we use the fact that a linear map is continuous and interpret open as relatively open. 



15 



any ball around Q is open. Thus, there is a point Q' in the interior which has a pre-image 
p' in the interior of T>{%^). We can find e G (0,1) and p" G 'Dj 1i j 2 ,...,j v ('Hn) such that 
p = (1 — e)p' + ep". Then the pre-images satisfy \\p — p'\\ = e\\p' — p"\\ < 2e. By choosing a 
sufficiently small ball around p one can make e arbitrarily small so that p, the pre-image of 
p, is in the interior of T>(Hn). 

To prove the second part, we need the fact that every interior point of P(^at) has full 
rank and consider the line p + pT with T G ker(T). Then T(p + pT) = T(p) = Q for all p. 
Since p has full rank, p + pT is an interior of T>{H^) for sufficiently small p. But since T 
is not positive semi-definite, the line p + pT must eventually cross the boundary of V{1-L N ) 
for some p > and some p < 0. 

In contrast to the anti-symmetric case, the symmetric subspace of %® N always has a non- 
trivial kernel for any choice of Q because all RDM for the pair of GHZ states |0 • • • 0)±| 1 ■ ■ ■ 1) 
are identical. 

Now recall that faces of V(T-Ln) correspond to subspaces of 1-Ln- One can apply Theorem 8 
with % N replaced by a subspace W N C Hn to obtain the following corollary. 

Corollary 9 If the map T : p h-> Q restricted to a face J 7 ofDil-L^) has a non-trivial kernel, 
then all interior points of T>j 1 j 2 ^ ^j v (J r ) have non-unique pre-images in J 7 . Moreover, every 
interior point ofT> Jx j 2 ^ j u [T^) has at least two pre-images on the boundary of T. 

When T maps two distinct pure states into the same Q, the line connecting them is also 
mapped to Q. Thus, if the RDM of a pure state do not determine it uniquely among all 
pure states, then there will always be mixed states with the same image. 

This leaves the following question. Suppose that \ip)(ip\ is the only pure state mapped 
to a RDM vector Q. Can there also be a mixed state p which is mapped to Q? By the 
discussion above, there must then be a line from |^)(^| to a state p' on the boundary which 
collapses to the single point Q under T. By assumption, p' can not be pure, but lies on 
a face J 7 of P("Hjv). Only in exceptional circumstances will T be one-to-one on T . Then 
by applying Corollary 9, one can find two additional mixed states p" and p'" which lie on 
the boundary of J 7 . Continuing this process inductively, we find that the pre-image of Q 
must be a cone whose apex is and for which typical extremal rays go from \i(j)(if)\ to 

a mixed state on a one-dimensional face (i.e., a line) on the boundary of T>(Hn). Thus, in 
general, if some mixed state has the same image as there will be a mixed state of 
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rank two with the same image. An equivalent characterization of rank two mixed states is 
that they are the reduction of pure state after extension to a qubit environment. 5 

In a series of papers, [8, 13, 14] Linden, et all showed that the answer to the question 
above is negative in certain special situations, including a triplet of qubits. Their main focus 
in the later papers is on the number of elements in R m needed to determine the state when 
the number of parties N is large. Their work uses parametric analysis so that "generic" can 
be interpreted as "almost every" in a probabilistic setting. Although settings with a trivial 
kernel seem extremely rare, we do not have a similar measure theoretic interpretation. 

V. SUMMARY AND OPEN QUESTIONS 

Let us summarize the various situations we know 

a) Almost all pure states are uniquely determined by a pair of RDMs pj, p'j which overlap 
in the sense that Jn J' ^ and JUJ' = {1,2,... N}. In the language of [5, 13, 14] we 
can say that almost all pure state are uniquely determined by a pair of RDM pjx, Pkl 
where J U K U L = {1, 2, . . . , N}; 

b) When both (pj and pjc) are extreme points of the set of RDM Uj(7-L N ) and Vjcil-L^) 
respectively, there is a unique pure state pre-image satisfying Trjc|^)(^| = pj and 
Trj|^)(^| = pjc; 

c) When 2m > N and p m = pi J i,...,m is an extreme point of the set of symmetric or 
anti-symmetric A-representable m-body RDM, it has a unique pre-image; 

d) When 2m < N, there are extreme points of the set of m-body consistent quantum 
marginals R m with non-unique pre-image. These map to extreme points of the set of 
m-body fermionic RDM with non-unique pre-image. For m = 2, the smallest N for 
which an example is known is N = 9. 

e) If the contraction map T : T>(T-Ln) ^ T>j 1 j 2 ^^j v {'Hn) has a non-trivial kernel, then all 
interior points of ^ , j 1 ,j 2 ,...,j„('Hv) have non-unique pre-images in V{1-L N ). Moreover, 

5 This fits into our framework using a subspace of %®( N + 1 ) for which Hn+i® C %® n (g> Q 2 , 



17 



in this situation, every interior point has at least two pre-images on the boundary of 
V{U N ). 

A number of open questions remain. 

A) Is there an extreme point of the set of A-representable m-body RDM or consistent 
m-body quantum marginals which is not exposed? 

B) What is the smallest N for which an extreme point of the set of 2-body consistent 
quantum marginals has a non-unique pre-image. Known results on QECC imply 
N > 5 [2]. 

C) If there is only one pure state with RDM vector Q, can there be a mixed state with 
the same RDM vector Q. We showed in the previous section that for generic situa- 
tions, any such mixed state has rank two or, equivalently, it suffices to consider qubit 
environments. 

The question of whether or not every face of the set of m-body reduced density matrices 
is exposed is of some theoretical interest. On the other hand, the relevance for the N- 
representability problem is not so clear. As pointed out by Kummer [12] the set of m-body 
A-representable RDM is the closed convex hull of its exposed points so that any extreme 
point which is not exposed is arbitrarily close to an exposed point. Moreover, for the original 
motivation of a constrained variational computation for an A-body Hamiltonian, it clearly 
suffices to consider the exposed points. 

Finally we remark that we use Hn Q %® n because that is the situation most commonly 
considered. However, most of our results apply more generally to subspaces of "H 1 ® H 2 ® 
■ ■ ■ <S> H N because such situations can be embedded as a subspace of some %® N . 

VI. ACKNOWLEDGEMENT 

JC is supported by NSERC and NSF of China (Grant No. 61179030). ZJ acknowledges 
support from NSERC and ARO. MBR is supported by US NSF. BZ is supported by NSERC 
and CIEAR. DZ is supported by NSF of China (Grant Nos. 10975181 and 11175247), and 



18 



NKBRSF of China (Grant No. 2012CB922104). 



[1] M. Altunbulak and A. Klyachko. The Pauli principle revisited. Comm. Math. Phys., 282:287- 
322, Sep. 2008. 

[2] A. Robert Calderbank, Eric M. Rains, P. W. Shor, and Neil J. A. Sloane. Quantum error 

correction via codes over GF(4). IEEE Trans. Inform. Theory, 44(4): 1369-1387, 1998. 
[3] M. Christandl and G. Mitchison. The spectra of quantum states and the Kronecker coefficients 

of the symmetric group. Comm. Math. Phys., 261:789-797, Feb. 2006. 
[4] A. J. Coleman. Structure of fermion density matrices. Rev. Mod. Phys., 35:668-686, Jul. 1963. 
[5] L. Diosi. Three-party pure quantum states are determined by two two-party reduced states. 

Phys. Rev. A, 70(1):010302, Jul. 2004. 
[6] R. M. Erdahl. The convex structure of the set of A^-representable reduced 2-Matrices. J. 

Math. Phys., 13:1608-1621, Oct. 1972. 
[7] C. Garrod and J. K. Percus. Reduction of the A^-particle variational problem. J. Math. Phys., 

5:1756-1776, Dec. 1964. 
[8] N. S. Jones and N. Linden. Parts of quantum states. Phys. Rev. A, 71:012324, Jan. 2005. 
[9] A. Klyachko. Quantum marginal problem and representations of the symmetric group. 

arXiv:quant-ph/0409113, Sep. 2004. 
[10] A. Klyachko. Quantum marginal problem and A r -representability. J. Phys.: Conf. Ser., 

36:72-86, Apr. 2006. 

[11] H. Kuhn. Linear inequalities and the pauli principle. Proc. Symp. Appl. Math., 10:141-147, 
1960. 

[12] H. Kummer. iV-representability problem for reduced density matrices. J. Math. Phys., 8:2063- 
2081, Oct. 1967. 

[13] N. Linden, S. Popescu, and W. K. Wootters. Almost every pure state of three qubits is com- 
pletely determined by its two-particle reduced density matrices. Phys. Rev. Lett., 89:207901, 
Oct. 2002. 

[14] N. Linden and W. K. Wootters. The parts determine the whole in a generic pure quantum 

state. Phys. Rev. Lett, 89:277906, Dec. 2002. 
[15] Y.-K. Liu, M. Christandl, and F. Verstraete. Quantum computational complexity of the 



19 



iV-representability problem: QMA complete. Phys. Rev. Lett., 98(11):110503, Mar. 2007. 

[16] S. A. Ocko, X. Chen, B. Zeng, B. Yoshida, Z. Ji, M. B. Ruskai, and I. L. Chuang. Quan- 
tum codes give counterexamples to the unique preimage conjecture of the A-representability 
problem. Phys. Rev. Lett, 106(11):110501, Mar. 2011. 

[17] M. Rosina. Some theorems on uniqueness and reconstruction of higher-order density matrices. 
In J. Cioslowski, editor, Many-electron densities and reduced density matrices. Springer, 2000. 

[18] M. B. Ruskai. iV-representability problem: particle-hole equivalence. J. Math. Phys., 
ll(ll):3218-3224, 1970. 

[19] M. B. Ruskai. Connecting iV-representability to Weyl's problem: the one-particle density 
matrix for N = 3 and R = 6 . J. Phys. A: Math. Gen., 40:F961-F967, Nov. 2007. 

[20] D. W. Smith. iV-representability problem for fermion density matrices. I. The second-order 
density matrix with A = 3. J. Chem. Phys. , 43:258, Jul. 1965. 

[21] C. N. Yang. Concept of off-diagonal long-range order and the quantum phases of liquid he 
and of superconductors. Rev. Mod. Phys., 34:694-704, Oct. 1962. 



